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ABSTRACT 

Munition  storage,  handling,  and  maintenance  facilities  consist  of  a  variety  of  structures. 
They  include  buried  rebar  re-enforced  concrete  walled  igloos  of  several  types,  as  well  as  above 
ground  buildings.  One  of  the  safety  hazards  of  concern  is  the  protection  of  these  munitions  from 
die  effects  of  lightning.  These  structures  are  electromagnetically  complex,  because  they  consist  of 
a  variety  of  inhomogeneous  materials  (e.g.,  concrete  with  rebar)  which  may  be  either  conducting 
or  partially  conducting.  In  addition,  the  structures  usually  have  metallic  penetrations  such  as 
electrical  cables  or  plumbing,  as  well  as  a  lightning  protection  system  including  an  earth  ground  of 
some  type. 

•  The  objective  of  this  paper  is  to  describe  how  the  lightning  hazards  to  such  structures  can 
be  evaluated  using  advanced  formulations  of  Maxwell’s  Equations.  The  method  described  is  the 
Three  Dimensional  Finite  Difference  Time  Domain  Solution.  It  can  be  used  to  solve  for  the 
lightning  interaction  with  such  structures  in  three  dimensions  and  include  a  considerable  amount  of 
detail. 


Examples  of  lightning  strikes  to  buried  igloos  and  above  ground  buildings  will  be 
presented.  The  physical  details  which  are  included  in  the  models  are  discussed.  The  results 
include  the  voltages  and  currents  induced  on  conductors  which  penetrate  the  facility,  as  well  and 
the  internal  electric  and  magnetic  fields.  Possibilities  for  internal  arcing  are  described.  These 
results  can  then  be  used  to  evaluate  the  possible  hazard  to  materials  stored  inside.  Of  special 
interest  is  the  evaluation  of  the  effectiveness  of  earth  ground  systems  and  how  they  affect  energy 
penetration  to  the  facility  interior. 

1.0  INTRODUCTION 

Lightning  hazards  to  explosives  storage,  handling,  and  maintenance  facilities  can  be 
described  in  terms  of  electric  and  magnetic  fields  (and  their  time  derivatives)  and  the  resulting  direct 
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and  induced  electrical  currents  which  are  found  in  and  around  critical  locations  of  the  facility  during 
a  lightning  strike. 

The  space  and  time  distributions  of  such  fields  and  currents  follow  solutions  of  Maxwells 
Equations  providing  that  appropriate  initial  and  boundary  conditions  can  be  supplied  in  the  regions 
of  interest  and  that  a  method  of  solution  can  be  applied. 

This  paper  describes  a  numerical  computer  model  which  applies  Maxwells  Equations  to 
describe  a  specified  lightning  attachment  to  a  specific  building  or  facility.  The  result  shows  how 
electromagnetic  fields  and  currents  are  distributed  in  space  and  time  in  and  near  the  facility  during 
the  lightning  strike. 

Examples  are  given;  1.  For  an  earth  covered  storage  igloo  with  iron  rebar  re-enforced 
concrete  walls,  and  2.  For  a  rectangular  building  with  cinder-block  walls  and  a  metal  roof.  Both 
structures  have  provisions  for  "lightning  protection"  in  the  form  of  air  terminals  connected  to  a 
ground  counterpoise  system.  It  will  be  shown  that  fields  and  currents  within  these  structures  can 
be  significantly  high  during  a  lightning  strike. 

2.0  DESCRIPTION  OF  THE  NUMERICAL  MODELS 


The  numerical  model  of  the  structure  and  surrounding  environment  is  based  upon  a  finite 
difference  time  domain  solution  of  Maxwell's  equations.  The  solution  technique  is  explicit  and 
accurate  to  second  order  in  the  time  and  spatial  increments,  which  in  these  models  correspond  to 
the  three  dimensional  cartesian  coordinate  increments  as  obtained  by  Merewether  and  Fisher  [1]. 

A  problem  space  containing  the  facility  and  surrounding  environment  is  divided  into 
rectangular  cells.  Each  cell  has  a  staggered  spatial  grid,  as  shown  in  Figure  1,  composed  of  the 
vector  components  of  E  and  H.  There  are  approximately  one  million  cells  in  the  lightning  strike 
problem  spaces  discussed  in  this  paper.  The  cell  dimensions  Ax,  Ay  and  Az  are  12”x6"x6"  for  the 
igloo  and  6"xl2"xl2"  for  the  building.  The  field  components  in  each  cell  are  calculated 
numerically  via  the  finite  difference  form  of  Maxwell’s  Equations  [1]. 


MAXWELL’S  EQUATIONS 


H®  +  VxE 

=  M 

(1) 

8^  +  oE-VxH 

=  -J 

(2) 

V  •  E 

=  P 
e 

(3) 

V  •  H 

=  0 

(4) 

Figure  1  Staggered  Spatial  Grid 
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The  time  step  (increment)  for  this  finite  difference  solution  of  Maxwell's  equations  is 
determined  by  the  Courant  criterion,  which  may  be  viewed  as  requiring  that  the  speed  of  numerical 
propagation  be  greater  than  the  fastest  physical  wave  speed,  in  this  case  the  speed  of  light  in  air. 
Specifically,  the  Courant  condition  is: 

At  <  — , - - - =-  (5) 

J  1  +  1  ♦  1 

V  Ax2  Ay2  Az2 

where  At  is  the  time  step.  Ax,  Ay,  and  Az  are  the  three  cartesian  spatial  increments  and  c  is  the 

speed  of  light  in  the  air.  For  the  igloo  At  is  .25  x  10’9  s  and  for  the  building  At  is  .33  x  10-9  s. 
The  smallest  spatial  increments  control  the  time  step,  but  the  largest  spatial  increments  determine 
the  bandwidth  of  the  solution.  The  rule  of  thumb  used  is  that  the  upper  frequency  limit  of  the 
solution,  fmax,  is  given  by 


fmax  -  •  (6) 

5  max  (Ax,Ay,Az) 

For  the  igloo  and  building  models  discussed  here,  this  corresponds  to  an  fmax  of  a  few  hundred 
MHz,  winch  is  more  than  sufficient  for  the  worst  case  lightning  environment  scenario. 

Maxwell's  curl  equations  (1),  (2)  form  a  system  of  hyperbolic  partial  differential  equations 
which  not  only  require  initial  conditions  at  all  spatial  locations,  but  also  the  boundary  values  of  the 
electromagnetic  field  components  (or  their  normal  derivatives)  at  all  times  to  obtain  a  well  posed 
solution.  These  values  must  be  supplied  at  the  boundaries  of  the  computational  volume  by  an 
appropriate  termination  condition.  The  boundary  condition  employed  was  derived  by  Mur  [2],  and 
is  essentially  a  first  order  integration  along  outgoing  (with  respect  to  the  interior  of  the 
computational  volume)  characteristics.  That  is,  the  characteristic  direction  is  chosen  to  be  causal  in 
time  and  along  the  outward  normal  to  the  bounding  surface,  which  is  a  two  dimensional  cartesian 
coordinate  plane.  Boundary  conditions  also  must  be  imposed  on  metallic  surfaces  such  as  the 
door,  interior  wall  and  metal  equipment  The  boundary  condition  on  metal  surfaces  at  least  as  large 
as  a  cell  face  is  that  the  tangential  electric  fields  at  the  surfaces  of  the  metallic  objects  are  set  equal  to 
zero  each  time  step.  Although  this  is  correct  only  for  perfect  electrical  conductors,  on  the  time 
scale  of  interest  it  is  an  excellent  approximation. 

If  the  Maxwell  divergence  equations  (3, 4)  are  satisfied  at  the  initial  time  step  then  the  finite 
difference  time  development  of  the  curl  equations  automatically  satisfy  the  divergence  equations  at 
each  time  step.  Thus  the  static  solution  in  the  problem  space  satisfying  (3)  and  (4)  is  tantamount  to 

specifying  the  initial  conditions  for  the  problem  The  simplest  initial  condition  is  to  set  E  =  H  =  p 
=  0  throughout  the  problem  space.  However,  physically,  a  lightning  discharge  is  normally  a 
dynamic  release  of  a  static  field  buildup  between  the  cloud  and  ground.  The  building  or  facility 
under  consideration  will  usually  cause  local  static  field  enhancements  from  the  charge  buildup 
between  cloud  and  ground.  The  air  dielectric  breakdown  will  then  usually  occur  at  the  point  of 
highest  electric  field,  e.g.,  an  air  terminal  or  protrusion  of  the  structure. 

Thus  it  is  sometimes  necessary  to  obtain  the  initial  static  solution  for  the  facility  under  high 
field  conditions  in  order  to  faithfully  track  the  fields  and  currents  of  the  resulting  lightning  strike. 
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At  other  times  it  may  be  sufficient  to  realize  that  under  linear  conditions  and  a  given  lightning 
current  injection  waveform,  the  final  solution  is  the  superposition  of  the  initial  static  solution  and  a 
dynamic  solution  with  the  initial  fields  and  charge  density  set  to  zero.  This  paper  will  be  primarily 
concerned  with  the  dynamic  part  of  the  solution  under  zero  initial  conditions. 

In  addition  to  the  appropriate  boundary  and  initial  conditions,  the  material  properties  at  each 
cell  location  must  be  specified.  This  consists  of  the  magnetic  permeability,  [X,  in  equation  (1);  the 

conductivity,  a,  in  equation  (2)  and  the  dielectric  constant,  e,  in  equations  (2)  and  (3).  If  the 
material  is  homogeneous  within  the  cell  (for  example,  volumes  of  air,  soil,  concrete,  etc.)  then  the 
appropriate  values  of  fi,  a,  and  e  are  included  in  the  time  advance  equations  for  the  cell  in  question. 

If  the  material  properties  are  inhomogeneous  in  each  cell  (detailed  structure,  etc.)  then  a 
decision  must  be  made  on  how  to  represent  the  properties  in  each  cell.  In  some  cases  average 
properties  are  sufficient  and  in  other  cases  they  are  not  Special  considerations  are  available  for 
treating  apertures  in  metal  walls  and  also  for  pipes  and  thin  wires  (radii  much  smaller  than  cell 
dimensions)  which  may  run  throughout  the  problem  space.  These  pipes  and  wires  can  be  carriers 
of  high  current. 

The  buildings  and  facilities  of  interest  usually  have  a  great  deal  of  "thin  wire"  situations  in 
the  form  of  signal  and  power  lines,  rebar  in  reinforced  concrete,  pipes,  plumbing,  metal  poles,  the 
lightning  protection  air  terminals,  down  conductors,  counterpoise,  etc. 

The  thin  wires  and  rods  are  implemented  in  a  self  consistent  fashion  by  making  use  of  the 
telegrapher's  transmission  line  equations.  The  telegrapher's  equations  (7),  (8)  are  a  one 
dimensional  solution  of  Maxwell's  in  terms  of  currents,  Iw,  and  voltages,Vw,  on  the  wires,  which 
are  required  to  have  diameters  less  than  cell  size  (spatial  increment).  The  per  unit  length 
inductances  and  capacitances  are  defined  (9),  (10)  with  respect  to  the  cell  size  and  the  wire 
diameter,  2a. 

One  Dimensional  Transmission  Line  Equations: 

^  =  -  Lw  -  IwRw  +  Ez(iwjwJc)  (7) 

Hz  =  "  Cw  ~3tL  ■  GwVw  (8) 

where  Lw  and  Cw  is  the  in-cell  inductance  and  capacitance  of  the  wire  per  unit  length. 


Cw  - 


27taeEr(a) 


'  w 


2jce 


(10) 


Gw  is  the  in-cell  conductance  from  the  wire  to  the  surrounding  conductive  medium 
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Gw  =  f  Cw  (11) 

The  wire  resistance  per  unit  length,  Rw,  is  obtained  by  considering  the  surface  conduction 
of  the  metal  in  question  using  the  skin  depth  obtained  for  a  frequency  of  1  MHz.  The  resistance 
for  pipes,  wire,  iron  rebar,  etc.,  is  normally  on  the  order  of  10"3  Ohms/meter.  In  practice,  the 
major  results  at  early  time  seem  to  be  relatively  insensitive  to  variations  of  the  resistance. 

In  the  computer  code,  the  wires  and  pipes  are  embedded  into  the  staggered  grid  and  are 
driven  by  the  electric  field  component  (see  equation  (7))  calculated  by  the  three  dimensional 
solution  of  Maxwell's  equations.  In  order  to  maintain  electrical  charge  conservation,  this  wire 
current  must  also  be  injected  back  into  the  driving  electric  field  component  as  a  source  current  via 
Maxwell's  Equation  (2).  At  the  interconnections,  which  are  voltage  nodes,  Kirchoff  s  law  is 
invoked.  At  locations  where  the  wires  are  situated  in  the  soil  or  concrete,  the  wires  are  in  electrical 
contact  with  the  soil  or  concrete  with  in-cell  conductance  given  by  Gw  in  equation  (11).  This  is 
also  true  of  the  facility  ground  wire  which  is  in  contact  with  the  soil. 

Complex  networks  of  thin  wires  (e.g.,  concertina  or  metal  rebar  mesh  embedded  in 
conducting  concrete)  are  included  in  the  model  by  a  vectorized  extension  of  the  transmission  line 
formalism.  Vectorized  average  wire  currents  coincide  with  the  electric  field  vectors  in  each  cell  and 
a  corresponding  average  inductance  and  resistance  is  associated  with  each  wire  current  vector.  Six 
component  tensors  exist  at  the  cell  comers  (nodes)  describing  the  equivalent  transmission  line 
voltages,  wire  capacitance,  and  conductance  to  the  embedding  medium.  A  36  component 
connectivity  tensor  exists  at  each  node  describing  the  ways  that  wires  are  connected  at  the  nodes. 

At  the  boundaries  of  the  problem  space,  some  termination  condition  must  be  applied  to  both 
the  counterpoise  extensions  and  the  power  and  signal  lines  and  metal  pipes  entering  the  problem 
space.  The  boundary  condition  is  applied  at  current  nodes  and  is  the  equivalent  of  the  Mur 
boundary  condition  applied  to  the  magnetic  fields  [2]. 

3,0  THE  LIGHTNING  STROKE  CURRENT  WAVEFORM  AND  INJECTION 

The  problem  is  initiated  by  imposing  a  pre-determined  lightning  wave  form  from  the  top 
edge  of  the  problem  space  to  a  specific  point  on  the  structure.  In  a  typical  computational  case 
described  below,  the  lightning  current  waveform  is  characteristic  of  a  1%  stroke  of  negative 
lightning.  The  lightning  current,  I(t),  is  given  as  a  function  of  time  by 


I(t)  =  1.1  x  105  sin2 
0  <_  t  <_  .5  x  10_6s 
I(t)  =  1.1  x  105  sin2 


ft  -  .5  x  10-6\ 


5  x  10 


-5 


A 


(12) 


.5  x  IQ’6  <  t 
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which  has  a  peak  current  of  110  kA  occurring  at  .5  fis.  The  lightning  current  appears  without 
propagation  delays  in  a  line  of  vertical  electric  fields  (Hz)  from  the  top  of  the  computational  volume 
to  the  attach  point.  The  lightning  current  is  injected  into  the  electric  fields  by  dividing  the  current 
by  the  cell  area  whose  normal  is  parallel  to  the  vertical  direction.  This  becomes  the  source  current 
density,  J,  in  Maxwell's  equation  (2).  A  number  of  different  parameters  are  studied:  lightning 
stroke  attachment  location,  soil  electrical  conductivity,  structure  wall  rebar  composition,  and  power 
box  attachment  at  the  walls  and  ceiling.  These  parameters  are  varied  in  order  to  provide 
environments  based  upon  the  range  of  situations  which  could  be  encountered. 

The  computer  model  contains  features  of  interest  such  as,  soil,  concrete,  rebar, 
counterpoise,  etc.,  which  are  included  in  the  computer  model  in  a  modular  form.  These  separate 
features  may  be  included  or  excluded  from  the  model  by  calling  subroutines  specific  to  the  features 
desired.  The  computations  are  performed  on  a  CRAY  II  computer.  Typical  run  times  are  1  hour 
of  computer  time  for  each  microsecond  of  real  time. 

4.0  LIGHTNING  STRIKE  MODELS 

The  analysis  of  the  preceding  sections  has  been  applied  to  two  structures:  (1)  an  earth 
covered  storage  igloo  with  iron  rebar  reinforced  concrete  walls  as  shown  in  Figure  2  and,  (2)  a 
rectangular  constructed  building  with  a  metal  roof  as  shown  in  Figure  3. 

The  igloo  interior  is  completely  surrounded  with  either  metal  or  iron  rebar  which  forms  a 
"leaky"  electromagnetic  shield  for  the  interior.  A  schematic  drawing  of  the  igloo  vertical  mid- 
cross-section  is  shown  in  Figure  4. 

The  building  is  made  of  concrete  block  outer  walls  with  no  rebar,  a  metal  roof,  and 
concrete  with  rebar  floor  and  inner  walls  with  rebar.  Thus  the  building  cannot  be  considered  as 
having  a  contiguous  shielding  effect. 

For  both  models  the  numerical  computer  output  from  a  simulated  lightning  strike  may  be 
categorized  as  follows: 

1 .  Contour  Plots  -  These  are  "snapshots  in  time"  of  the  electric  and  magnetic  field 
structures  on  a  plane  cross-section  of  the  building  at  some  time  after  the  initiation  of 
the  strike. 

2.  Time  Dependent  Plots  -  these  are  time  dependent  graphs  of  electric  and  magnetic 
fields  at  selected  points  in  the  problem  space.  Currents  and  voltages  on  thin  wires 
and  rods  also  have  time  dependent  plots  at  selected  points. 

3 .  Current  Arrays  -  These  are  spreadsheet  tabulations  of  wire  currents  in  specific  areas 
of  the  building. 

4.  Field  Maxima  -  These  are  computer  searches  at  selected  times  to  find  the  maximum 
electric  and  magnetic  fields  and  the  maximum  time  derivative  of  the  magnetic  field 
within  a  specified  boundary  inside  the  building. 

Figure  5  shows  a  contour  plot  of  the  vertical  mid-plane  longitudinal  cross-section  of  the 
igloo  corresponding  to  the  schematic  in  Figure  4.  The  electric  field  pattern  outlines  some  of  the 
prominent  features  of  the  igloo,  i.e.,  the  z-cage,  soil  berm  over  the  igloo,  headwall,  backwaU,  etc. 
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Figure  3  Building  -  Right  Side  View 
Protection  System 


;loo  --  Lightning  Strike  Model 


Window  Screens  and  Lightning 


Figure  4  Igloo  Vertical  Cross-Section  at  j  =  jm  =  75 


Figure  5  Electric  Field  Vector  and  Magnitude  Contour  Plot  for  Vertical  Mid- 
Cross-Section  of  Igloo 
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The  vectors  show  the  projection  of  the  electric  field  vector  at  each  cell  onto  the  mid-plane  at  a  time 
1  |isec  after  the  initiation  of  the  strike.  The  length  of  the  vector  is  proportional  to  the  logarithm  of 
the  electric  field.  The  contour  lines  show  lines  of  equal  electric  field  magnitude  labeled  as  powers 
of  10  of  the  field  magnitude  in  volts  per  meter.  For  example,  the  line  labeled  4.0  represents  field 
magnitudes  of  10,000  volts/meter. 

Figure  6  shows  a  contour  plot  on  a  vertical  x-z  plane  of  the  building  cutting  through  wire 
mesh  on  the  window  nearest  the  strike.  The  view  is  as  if  looking  from  the  back  of  the  building. 
The  field  patterns  show  essential  geometrical  features  of  the  model,  i.e.,  roof,  supporting  I-beams, 
outer  wall,  etc. 

The  window  mesh,  a  wire  grid  covering  the  building  windows,  is  being  charged  (note  E- 
field  vectors  pointing  away  from  the  mesh)  and  appears  to  focus  the  electric  field  into  the  interior  of 
the  building.  The  field  levels  are  very  high  within  the  building  approaching  1  Megavolt/meter 
(contours  are  labeled  as  powers  of  10  of  the  electric  field  magnitude). 

In  this  case.  Figure  6,  the  lightning  protection  system  is  not  connected  to  the  metal  roof.  At 
.462  p.  seconds  the  top  of  the  roof  is  positively  charged  and  the  bottom  of  the  roof  is  negatively 
charged. 

Figure  7  shows  the  effect  of  adding  an  I-beam  (perpendicular  to  the  contour  plane)  with  a 
hanging  metal  cable  hoist.  The  field  at  the  bottom  of  the  hoist  is  on  the  order  of  a  few 
megavolts/meter  and  represents  a  potential  for  arcing  between  the  hoist  and  the  floor  rebar  (or  any 
other  piece  of  grounded  equipment).  In  this  case  the  lightning  protection  system  is  in  contact  with 
the  metal  roof  which  is  also  in  contact  with  the  I-beam. 
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Figure  8  shows  time  dependent  plots,  corresponding  to  Figure  7,  of  the  lightning  injection 
current  (given  by  equations  (12)),  the  electric  field  and  wire  voltage  in  the  middle  of  the  window 
screen,  and  the  voltage  between  the  hoist  hook  and  the  floor  rebar.  This  is  a  case  showing  that 
connecting  the  lightning  protection  system  to  the  building  structure  can  enhance  the  hazard  inside 
the  building. 

5.0  INDUCTIVE  AND  CAPACITIVE  COUPLING  TO  THE  INTERIOR  OF  THE 

IGLOO 

The  construction  of  the  igloo  provides  that  the  interior  of  the  igloo  is  completely 
surrounded  by  a  "leaky"  electromagnetic  shield  consisting  of  rebar  in  conductive  concrete  and 
metal  doors,  walls,  etc.  If  is  of  interest  to  examine  the  character  of  electromagnetic  energy  leaking 
into  the  igloo  interior  from  the  point  of  view  of  the  model. 

Figure  9  shows  time  dependent  plots  the  lightning  injection  current  waveform  at  the  air 
terminal  and  electric  field  components  at  a  point  on  the  igloo  center  line  near  the  back  wall  and 
ceiling  of  the  igloo.  In  this  case  the  igloo  contains  only  the  internal  metal  wall.  The  strong  Ex 
component  of  the  field  peaks  at  .5  fisec  and  its  waveform  appears  to  follow  the  time  derivative  of 
the  lightning  current  injection  waveform.  This  is  interpreted  to  mean  that  there  is  an  inductive 
coupling  between  currents  flowing  on  longitudinal  elements  of  rebar  and  the  interior 
electromagnetic  field. 
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Figure  6 


Electric  Field  Vector  and  Magnitude  Contour  Plot  for  a  Vertical  Plane 
Passing  Through  the  Window  Mesh  of  the  Building 
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Figure  7 


Hectric  Field  Vector  and  Magnitude  Plot  for  Building  Showing  the 
Effect  of  an  Internal  I-Beam  and  Metal  Cable  Hoist 


862 


VOLTAGE  kv  CURRENT  Kfi 


ey  kv/w  CURRENT  KO 


13 


In  contrast,  the  vertical  Ez  field  peaks  at  1  nsec  and  its  waveform  is  similar  to  the  lightning 
current  injection  waveform  except  that  it  decays  more  rapidly.  This  is  presumably  capacitive 
coupling  due  to  charge  collecting  on  the  rebar  which  is  in  contact  with  the  lighting  protection 
system  at  several  different  grounding  locations.  The  decay  in  field  strength  at  late  time  appears  to 
be  due  to  charge  leaking  off  the  rebar  onto  the  counterpoise  system  and  by  conduction  through  the 
concrete  and  soil. 

The  inductive  and  capacitive  coupling  is  illustrated  globally  by  comparing  interior  field 
contour  plots  at  .5  jisec  (Figure  10)  and  1  (isec  (Figure  1 1).  In  Figure  10  the  Ex  component  is 
large  throughout  much  of  the  interior  and  at  later  time  (Figure  1 1)  the  vertical  Ez  field  dominates. 

The  charge  collecting  on  the  rebar  may  be  noted  in  Figures  10  and  1 1  by  observing  electric 
field  vectors  pointing  away  from  the  ceiling  and  back  wall  in  both  directions. 

It  is  noted  that  the  largest  fields  are  near  the  floor  in  both  figures  and  are  on  the  order  of 
100  Kvolts/meter  at  1  |isec.  This  is  due  primarily  to  capacitive  coupling  of  charge  on  the  rebar 
which,  again,  is  enhanced  by  electrical  contact  between  the  lightning  protection  system  and  the 
igloo  metal  structure. 

6.0  CONCLUSIONS 

A  numerical  computer  model  of  Maxwell's  Equations  has  been  applied  to  buildings  typical 
of  munitions  storage  and  handling  structures  to  calculate  potential  hazards  due  to  lightning  stakes. 
It  is  seen  that  detailed  electromagnetic  field  profiles  and  currents  may  be  calculated  which  estimate 
in  a  realistic  manner  the  hazardous  areas  in  and  around  the  facility.  The  possibility  of  both 
inductive  and  capacitive  electromagnetic  coupling  to  the  interior  of  the  structures  has  been 
demonstrated.  This  coupling  can  be  enhanced  by  electrical  contact  between  the  lightning  protection 
system  and  the  metal  structural  components  of  the  facility.  Hazardous  non-linear  effects  such  as 
electrical  arcing  and  explosive  decomposition  of  building  structure  (e.g.,  shrapnel  from  pieces  of 
exposed  wood  or  concrete)  are  possible  for  the  calculated  examples.  These  results  and  techniques 
may  be  applied  to  evaluate  potentially  hazardous  explosive  storage  and  handling  situations. 
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